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In this paper, the dynamical heterochnic orbit and attractor have been employed to make the 
late-time behaviors of the model insensitive to the initial condition and thus alleviates the fine 
tuning problem in cosmological dynamical system of barotropic fiuid and quintessence with a double 
exponential potential. The late-time asymptotic behavior of the double exponential model does 
1^-^ \ not always correspond to the single case. The heteroclinic orbits are the non-stationary solutions 

. and in general they will interpolate between the critical points. Although they can not be shown 

' analytically, yet a numerical calculation can reveal most of their properties. Varied heteroclinic 

I orbits and attractors including tracking attractor and de Sitter attractor have been found. 
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I. INTRODUCTION 

> , - n 

Q^ \ Astronomical observation on the cosmic microvifave background(CMB) anisotropy IJ, supernova type Ia(SNIa)|3| 

and SLOAN Digital Sky Survey(SDSS) 3] depicted that our Universe is spatially flat, with about seventy percent of 
■ the total density resulting from dark energy that has an equation of state w < —1/3 and accelerates the expansion 
of the Universe. Several candidates to represent dark energy have been suggested and confronted with observation: 
cosmological constant, quintessence with a sing le field3 or with N coupled field 's'l, phantom field with canonical f6j 
^5 , or Born-Infield type Lagrangianj2|, k-essence,Si] and generalized Chaplygin gas(GCG) 9]. 

I One of the most important issues for dark energy models is the fine tuning problem, and a good model should 
limit the fine tuning as much as possible. The dynamical attractor of the cosmological system has been employed to 
make the late time behaviors of the model insensitive to the initial condition of the field and thus alleviates the fine 
O tuning problem. In quintessence |0| and phantom fTJ models, the dynamical systems have tracking attractors that 
^ , make the quintessence and phantom evolve by tracking the equation of state of the background cosmological fluid so 
^ ' as to alleviating the fine tuning problem. In addition, there are also two late time attractors in the phantom system 
. . corresponding to the big rip phase 12] and de Sitter phase'l3|. On the other hand, exponential potentials can arise 
from string/A/-theory, e.g. via compactification on product spaces possibly with fluxes. In this case, the equations 
of motion can be written as an autonomous system, and the power-law and de Sitter solutions can be determined 
^ by an algebraic method. The properties of the attractor solutions of exponential potentials can lead to models of 
' quintessence jl^. And general exact solution for double exponential potential with one exponent is the negative of 
the other for quintessence was presented in Ref.[l^. The aim of this paper is to study the cosmological dynamics 
of barotropic fluid and scalar field with a double exponential potential and point out that the late-time asymptotic 
behavior does not always corresponding to the single-exponential caseJii]. We show that the existence of tracking 
attractor and de Sitter attractor. We also flnd some heteroclinic orbits which mean solutions interpolate between 
different critical points. Emphasis must be placed on that the tracking orbits have the similar but not exactly equal 
dynamical behavior and the initial possibilities consist in a wide range. 

II. PHASE SPACE AND CRITICAL POINTS 

We consider 4-dimensional gravity with barotropic fluid and a scalar (j) which only depend on cosmic time t. The 
scalar has a double exponential potential 

-ai(p 



U(^) = Aie-^i-^ + Aae-"^-^. (1) 
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Since current observations favor fiat Universe, we will work in the spatially flat Robertson- Walker metric. The 
corresponding equations of motion and Einstein equations could be written as 

^2 

H = — -{p-^ +P7 + 0^), 
Py = -3H{p^+p^), 
4> + 3H^ + V'{(f>) = 0, 

i^' = y(P7 + P0), (2) 

where — SttG, p-^ is the density of fluid with a barotropic equation of state = (7 — l)p^, where < 7 < 2 
is a constant that relates to the equation of state by w = 7 — 1. The overdot represents a derivative with respect 
to t, the prime denotes a derivative with respect to (p. p^ — ^cfP + V{(j)) and — — V{(j)) are the energy 
density and pressure of the (j) field respectively, and H is the Hubble parameter. Phase space methods are particularly 
useful when the equations of motion are hard to solve analytically for the presence of barotropic density. In fact, the 
numerical solutions with random initial conditions are not a satisfying alternative because of these may not reveal 
all the important properties. Therefore, combining the information from the critical points analysis with numerical 
solutions, one is able to give the complete classification of solutions according to their late-time behavior. Similarly 

as in Ref.fiol fl^. we introduce the following dimensionless variables x — U — '"^^^ ^ ~ '"^^h 

— ^'•'^•'^ ("i^) a,nd N = logo. Then, the Eqs.© could be reexpressed as the following system of equations: 



z[7(l -x^ -y^ - z^) + 2x^\ - -\\'^aixz. (3) 



Also, we have a constraint equation 

+ (4) 

where 

Different from the case in the single exponential potential, the parameter F is dependent on 

(aij/^ -I- 0122:^)^ 

The equation of state for the scalar field could be expressed in terms of the new variables as 



(6) 



P0 x'^ -y^ - z'^ 

w^^ — = ^- — — — J. (7) 

Pcj> x^ + z^ 

According to the Eqs.lO, one can obtain the critical points and study the stability of these points. Substituting 
linear perturbation x = x-\- 6x^ y = y + Sy, z = z + Sz near the critical points into the three independent equations, to 
first order in the perturbations, gives the evolution equations of the linear perturbations, from which we could yield 
three eigenvalues. Stability requires the real part of all eigenvalues to be negative. The results are contained in Table 
I. 
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Case Critical points (x,y,z) 







Stability 


a'f 


7-1 


stable 


"2 


-y - 1 


stable 


1 


2 


stable 


1 




stable 


1 


-1 


stable 


1 


1 


unstable 



V2~'V2V ' 

W \/2^'U'V2\/ ^1 



1/2 /-I _ a2>,-l/2 



(v) 0,(l-fJ)-^/^(l-f^ 

(vi) ±1,0,0 



TABLE I: The properties of the critical points. 



Cases (i) and (ii): These are tracking attractors. The hnearization of system © about these critical points yields 
three eigenvalues: 



■3(ai-Q;2)7 3{Q!f(7 
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-2) 


- 2472^2]} 
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- V«i(7-2)[a2(97 


-2) 


- 2472^2]} 






4a2 






3{«i(7- 


2)4 


V«i(7-2)K(97- 


2)- 


2472^2]}^ 



4a2 



(8) 



(9) 



According to the eigenvalue (8), the stability of the attractor requires the condition ^/Stk < ai < a2 or 02 < ai < 
—y/3jK. For the eigenvalue (9), it requires the condition y/3jK < a2 < ai or ai < a2 < — v^^tk. 

Cases (iii) and (iv): These are quintessence attractors. The linearization of system about these critical points 
yields three eigenvalues: 



ai{ai — a2) al 
2^2 

a2{ct2 — Oil) al 



37 



2k2 



(10) 
(11) 



According to the eigenvalue (10), the stability of the attractor requires the condition < ai < \/3^k and ai < a2 
or —y/SjK < ai < and ai > a2- For the eigenvalue (11), it requires the condition < 0:2 < a/Stk and a2 < ai or 
— y/3jK < q;2 < and a2 > ai. 

Cases (v): The critical point is a attractor corresponding to equation of state w = —1 and cosmic energy density 
parameter — 1, which is a de Sitter attractor. The condition of such a de Sitter attractor is aia2 < 0. The 
linearization of system ||2J) about these critical points yields three eigenvalues: 



-37, 



-3^2 - Vl2aia2K2 +9k4 -3^2 + Vl2aia2K2 + 9k'' 



2^2 



2k2 



(12) 



Cases (vi): These critical points corresponding to kinetic-dominated solutions in the asymptotic regime. The 
linearization of system Q about these critical points yields three eigenvalues: 



(13) 
(14) 



Although the kinetic-dominated critical points are always unstable, we will find some new properties in numerical 
analysis. 
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FIG. 1: The attractor property of the quintessence at the presence of dust matter. We choose ai = 2, 02 = \/5, 7 = 1, 
K = 1. The heterochnic orbit connects the critical point which corresponds to the case (i) to the kinetic-dominated critical 
point (1,0,0). 





FIG. 2: The attractor property of the quintessence at the presence of dust matter. We choose qi = 1, 02 = 2, 7 = 1, k = 1. 
The heterochnic orbit connects the critical point which corresponds to the case (iii) to the kinetic-dominated critical point 
(1,0,0). 



III. HETEROCLINIC ORBITS 



Critical points are always exact constant solution in the context of autonomous dynamical systems. These points 
are often the extreme points of the orbits and therefore describe the asymptotic behavior. If the solutions interpolate 
between critical points they can be divided into heterochnic orbit and homoclinic orbit (closed loop). The heterochnic 
orbit connects two different critical points and the homoclinic orbit is an orbit connecting a critical point to itself. 
If the numerical calculation is associated with the critical point analysis, then we will find all kinds of heterochnic 
orbit, as show in Figure. 1 to Figure. 3. The initial point is (±1, 0, 0) for all orbits since it is a repeller. Especially, the 
heterochnic orbit is shown in Figure 2, which connects the tracking attractor to (1, 0, 0). In Figure 3, the heterochnic 
orbit connects the de Sitter attractor to (1, 0, 0). Tracking behavior consists in the possibility to avoid the fine tuning 
of initial conditions, obtaining the similar behavior for a while range of initial possibilities. If we take cosmic time 
t = to, the phase space contain a two-dimensional submanifold corresponding to set of initial possibilities. This initial 
submanifold consists with the intersection set of heterochnic orbit and t = to surface. 

In Figure 4, we plot the dynamical evolution of matter, radiation and dark energy for the model with double 
exponential potential, in which tracker mechanism was used to provide a non-negligible energy density at early 
eD0ch|l7j. The dynamical evolution of equation of state w va N are given in Figure 5. 



FIG. 3: The attractor property of the quintessence at the presence of dust matter. We choose ai = —1, Q2 = 1, 7 = 1, 
K = 1. The heterocUnic orbit connects the critical point which corresponds to the case (v) to the kinetic-dominated critical 
point (1,0,0). 




FIG. 4: Evolution of cosmic parameters for matter, radiation and quintessence. Plot begins from the equipartition epoch. 




FIG. 5: The evolution of the equation of state w vs N. We choose ai = —2, —1.5, —1, —0.5 respectively and a2 = 1, 7 = 1, 

K = l. 
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IV. CONCLUSION AND DISCUSSION 

In this paper, we investigate the cosmological dynamics of scalar field with a double exponential potential. We show 
that the late-time asymptotic behavior does not always corresponding to the single-exponential case. Furthermore, 
we study varied heteroclinic orbits including the tracking one. Tracking behavior consists in the possibility to avoid 
the fine tuning of initial conditions. The numerical calculation show that the tracking orbits have the similar but not 
exactly equal dynamical behavior and the initial possibilities consist in a wide range. 

Why there are new features for double exponential potential? First of all, the parameter F here is dependent on the 
scalar field in contrast with single case. The second reason is that new features will occur in the regime where both 
exponential terms are remarkable. The third possible reason is that the heteroclinic orbits connect not only late-time 
but also early-time behavior. 

We can extend our discussion to the O(A^) quintessence model0, the initial energy density could be larger and would 
not affect the evolution too much because the 'angular' kinetic energy decreases with a"''. In a scalar field model, 
tracker mechanism was used to provide a non-negligible energy density at early epoch. In the 0{N) quintessence 
model, if equipped with the same tracker mechanism, it will admit a wider range of initial energy density 
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